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Abstract 

The multiobjective design of digital filters using spiral optimization technique is considered in this paper. This 
new optimization tool is a metaheuristic technique inspired by the dynamics of spirals. It is characterized by its 
robustness, immunity to local optima trapping, relative fast convergence and ease of implementation. The 
objectives of filter design include matching some desired frequency response while having minimum linear 
phase; hence, reducing the time response. The results demonstrate that the proposed problem solving approach 
blended with the use of the spiral optimization technique produced filters which fulfill the desired characteristics 
and are of practical use. 

Keywords: Multiobjective filter design; Spiral optimization technique; Magnitude response; Minimum linear phase; 
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Introduction 

Digital filters exist in two types: Finite impulse response 
(FIR) and Infinite impulse response (IIR) or recursive. 
FIR filters suffer from the problem of high order (hence 
implementation and performance issues) if strict require- 
ments are imposed at the design stage. Furthermore, IIR 
filters can have smaller group delay than its equivalent 
FIR filters (Antoniou 1993, Antoniou 2005). The optimal 
design of an infinite impulse response (IIR) filter con- 
sists in choosing a set of coefficients of the filter to have 
a frequency response that optimally approximates the de- 
sired response (Antoniou 2005; Dumitrescu and Niemisto 
2004; Lai 2009; Tseng and Lee 2002, Tseng 2004; Ho 
et al. 2008; Quelhas and Petraglia 2009; Sanathanan and 
Koerner 1963). 

Different techniques exist for the design of digital fil- 
ters. Windowing method; in which the ideal impulse re- 
sponse is multiplied by a window function, is the most 
popular. There are various kinds of window functions 
(Butterworth, Chebyshev, Kaiser etc.), depending on the 
requirements on ripples in the passband and stopband, 
stopband attenuation and the transition width. These 
various windows limit the infinite length impulse response 
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of ideal filter into a finite window to design an actual 
response. Furthermore, windowing methods do not allow 
sufficient control of the frequency response in the vari- 
ous frequency bands and other filter parameters such 
as transition width. The designer always has to com- 
promise between the design specifications (Antoniou 1993, 
Antoniou 2005). 

Due to the presence of the denominator of the transfer 
function, the stability condition of the filter should be 
taken into account in the optimal design (Antoniou 1993, 
Antoniou 2005; Lu et al. 1998; Saab et al. 1999; Tseng and 
Lee 2002; Omoifo and Hinamoto 2004; Tseng 2004; Jiang 
and Kwan 2009, Jiang and Kwan 2010a; Lai and Lin 2010), 
resulting in a constrained optimization problem. Several 
sufficient conditions (Lang 2000; Lu 2000; Antoniou 2005; 
Ho et al. 2008; Pan 2009; Jiang and Kwan 2010b) have been 
established for the parameterization that represents the fil- 
ter s denominator by a single polynomial. The triangle- 
based stability conditions (Antoniou 1993) are necessary 
and sufficient and have been incorporated into several de- 
sign procedures (Lu 1998; Lu and Hinamoto 2003; Lu 
2006) that formularize the filters denominator by cas- 
caded second-order sections (SOSs). In (Lu 1998), vari- 
able transformation is used to convert the finite stability 
region into the entire coefficient space, such that the 
original constrained design problem becomes an un- 
constrained one in the transformed space. However, the 
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transformation increases the nonlinearity of the ob- 
jective function, which makes it hard to find good (global 
optimum) solutions in general In (Lu and Hinamoto 
2003), a perturbed stability triangle is proposed to guar- 
antee the SOS to have its zeroes inside a circle of given 
radius. It is combined with the Gauss-Newton strategy, 
resulting in an improved design. In (Lang 1999), the con- 
ditions presented for the SOS with zeros inside a circle of 
given radius enclose a triangular stability domain and can 
be easily incorporated into any constrained opti- 
mization formulations based on the SOS paramete- 
rization. A method that divides the overall design of an IIR 
filter into successive designs of its second-order sections is 
presented in (Saab et al. 1999), where one section is first 
designed, and then, another section is appended until all 
sections are designed. 

Because of finite word length effects occurring in prac- 
tical implementations of the designed filters, not only 
stability of the filter is of great importance but a stability 
margin is necessary as well. The poles of the transfer 
function should not lie too close to the unit circle. The 
sensitivity of pole locations to coefficient quantization 
increases with decreasing distance from the unit circle. 
Poles close to the unit circle may considerably enhance 
quantization noise and increase the maximum amplitude 
of small scale limit cycle. Consequently, it is desirable to 
have control over the maximum radius when designing 
IIR filters. 

Linear-phase filters are usually designed as non-recursive 
(FIR) filters which can have constant group delay over the 
entire base-band. However, when highly selective filters are 
required, a very high filter order is needed which makes 
these filters uneconomical or impractical. To eliminate this 
problem, attempts have been made to develop methods to 
design recursive (IIR) filters whose delay characteristics ap- 
proximate a constant value in the passband. This includes 
IIR filter design approach that can satisfy both magnitude 
and phase characteristics simultaneously (Inukai 1980; 
Cortelazzo and Lightne 1984; Sullivan and Adams 1997; 
Lang 2000; Lertniphonphun and McClellan 2001). The 
design of IIR filters with constant group delay in the 
passband is also carried out by using allpass structures 
through evaluation of phase response instead of ap- 
proximating the group delay directly (Jing 1987; Ikehara 
et al. 1992; Lang and Laakso 1994; Zhang and Iwakura 
1999). Some other methods used an indirect approach 
based on model reduction techniques where a linear- 
phase FIR filter that meets the required specifications 
is first designed and then a lower order IIR filter that 
meets the original amplitude specifications while main- 
taining a linear-phase response in the passband is obtained 
(Sreeram and Agathoklis 1992; Peng et al. 1992; Beliczynski 
et al. 1992). Coretlazzo and Lightner (Coretlazzo and 
Lightner 1984) have achieved the simultaneous design in 



both magnitude and group delay of IIR and FIR fil- 
ters based on multiple criterion optimizations. Lutova 
(Lutova 1997) has developed a new design method for el- 
liptic IIR filters that provide the implementation of half of 
the multiplication constants with few shifters and adders. 
Sullivan et al. (Sullivan James and Adams 1998) have pro- 
posed the algorithm based on the peak-constrained least- 
squares optimality criterion for cascaded IIR filters, which 
can design a filter that has an equalized group delay with- 
out the use of all pass filters, and it can simultaneously 
meet the frequency response magnitude specifications by 
using all of the filter coefficients available to optimize the 
filter. Lang (Lang 2000) has used least square method for 
designing IIR filter with prescribed magnitude and phase 
response. This parameterization of the transfer function 
has been used for designing IIR filters. Gordanajavanovic 
(Gordana 2006) has proposed a method for the design 
of IIR notch filters with desired magnitude characteris- 
tic, which can be either maximally flat or equi-ripple. Xi 
Zhang (Zhang 2008) have proposed a novel method 
for designing maximally flat IIR filters with flat group delay 
responses in the pass-band. 

Under these circumstances, evolutionary and metaheu- 
ristic optimization methods find their place. These are re- 
ferred to as global optimizers while the more familiar, 
traditional techniques such as conjugate gradient and the 
quasi-Newtonian methods are classified as local opti- 
mizers. The distinction between local and global search of 
optimization techniques is that the local techniques pro- 
duce results that are highly dependent on the starting point 
or initial guess, while the global methods are totally inde- 
pendent of the initial conditions (Recioui 2012). Though 
they possess the characteristic of being fast in convergence, 
local techniques, in particular the quasi-Newtonian tech- 
niques have direct dependence on the existence of at least 
the first derivative. In addition, they place constraints on 
the solution space such as differentiability and continuity; 
conditions that are hard or even impossible to satisfy in 
some situations (Recioui 2012). 

Previously, global optimization techniques have been 
implemented in the design of digital filters. One such 
approach using neural networks has been described in 
(Wang et al. 2006). Also, use of PSO in the design of 
frequency sampling finite impulse response (FIR) filter 
has been described in (Wang et al. 2004; Krusienski and 
Jenkins 2004; Chen and Chen 2006). Differential evolu- 
tion has been used in the design of digital filters has 
been implemented in (Karaboga D 2004; Storn 1996, 
2005; Karaboga N 2005). (Kit Sang and Kim-Fung 1998; 
Karaboga N 2004) have used Hierarchical Genetic Algo- 
rithms to the design and optimization of IIR filter struc- 
tures. Use of Particle Swarm Optimization (PSO) and 
Genetic Algorithms (GA) in the design of digital filters is 
described in (Ababneh, and Bataineh 2007). 
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In this work, the application of the novel optimization 
technique called spiral optimization to the design of 
digital filters is considered. The purpose is to design a 
filter that can simultaneously satisfy multiobjective cri- 
teria including frequency response and linear phase with 
the least possible group delay. 

Problem formulation 

Digital filters find their applications in different areas. 
One area is power system protection where measure- 
ment systems involve faulted signals associated with DC 
decaying signals, harmonic and sub-harmonic compo- 
nents. To eliminate these unwanted components, a digital 
filter design based on multi-objective optimization tech- 
nique to satisfy different specifications such as high speed 
response for a real-time application and frequency domain 
requirements. 

Digital filtering approach 

A digital filter based solution is proposed to remove un- 
wanted disturbances using digital filter design tech- 
niques. The filter time response must be included in the 
requirements. The present filtering application imposes 
different kind of specifications. On one hand, the time 
domain requirement where both a high speed and accur- 
ate system response are needed. On the other hand, the 
frequency domain requirements (DC, sub-synchronous 
and harmonic components elimination) which are the 
magnitude response within small bandwidth including 
sharp frequency edges as well as an approximately con- 
stant group delay in this band are required too. Usually 
the best optimum value of all the objective functions of 
this filter design can be obtained for some values of de- 
sign variables. A compromise or a trade-off between the 
objective functions must be made to achieve a satisfac- 
tory filter design. 

The considered recursive digital filter must satisfy 
three multi-objective functions. These functions are: 
1) meet a specified or a desired magnitude response spe- 
cification; 2) an approximately constant group delay; and 
3) a minimum time response or settling time which in- 
volves a minimum phase or a group delay. The optimi- 
zation approach considers the discrete-time transfer 
function which is formulated on the basis of some desired 
amplitude response and a stability margin parameter. A 
norm of the weighted error function is then minimized 
with respect to the transfer-function coefficients with a 
prescribed maximum pole radius referred to as stability 
margin. The stability margin parameter is varied to 
optimize the filter coefficients which minimizes mainly the 
magnitude response, satisfies the best approximately con- 
stant group delay and the lowest group delay that leads to 
minimum settling time or time delay of the system dy- 
namic response. 



Filter transfer functions 

In the general case an IIR filter can be described by its 
discrete-time difference equation 



y[n] + c iy[ n ~ l ] = Y d A n ~j] 

i=l /=0 



(i) 



Where x[n] and y[n] are discrete-time input and out- 
put signals. Equation (1) can be transformed into the 
Z-domain and assuming q and d[ are real coefficients a 
second order form transfer function can be obtained, 
having 2 M conjugate zeros and 2 N conjugate poles; 
called second order sections (SOS), as: 



H {z) =H 0 n ^ + ^ + Z l 

W U^bof + byZ + Z 2 



(2) 



Where a u and bjj are real coefficients and H 0 is a posi- 
tive multiplier constant. The polar formulation is also 
useful and is written as: 



H(z) = H( 



n™(z-r af <j e <)(z-r a) e-*><) 
n^{z-r bi ^>){z-r bi e-^) 



(3) 



Where r a j, 9 a j and r b j,6bj are the radii and angles of the 
zeros and poles, respectively. 

Filter stability margin 

As the poles are moved toward the origin (decreasing 
the pole radius), the system stability margin param- 
eter increases and the system settling time decreases. 
This action, in fact, brings two required and import- 
ant properties to the designed system. First, the sys- 
tem time or dynamic response is enhanced as settling 
time is decreased. Second, the system stability becomes 
more robust which is a very useful property, particu- 
larly in practical implementation. Indeed, the rounding 
or truncation of the filter coefficients may lead to an 
unstable implementation if the stability margin is too 
small. It is therefore desirable to approximate a given 
response by a transfer function with a prescribed max- 
imum pole radius named stability margin as shown in 
Figure 1. 

Magnitude response objective function 

The amplitude and the phase responses of a recursive 
filter is given by 

M(x, co) = \H(<j»% <p(x, co) = arg(tf (4) 
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between the actual amplitude response M(x, o) and the 
Where co is the frequency and x is a column vector desired amplitude or magnitude response M d (o>) as: 



with 2 M + 2 N + 1 elements and depend on the used 
formulation, that is a Cartesian form as 



x = [aoia n a2iboib n b2i...b2MHo] 
Or in polar form as 

* = [paPb H o] 

Where 



Pa 



And 



Vol, "ali-", r aMi "aM\ 



Pb = [ r bl, Obi, rbN: ObAl} 



(5) 



(6) 



(7) 



(8) 



The superscript T denotes the transpose operation. An 
approximation error can be formulated as the difference 



e(x) = M(x,(o)-M d (a>) 



(9) 



By sampling the error function e(x,o)i), the actual and 
the desired amplitude responses M(x,cOi) and Ma(co) at 
frequencies: coi, co 2 , co^, the column error vector is 
constructed as: 



[£(#)] = [ei(x),e 2 (x), -^e^x)] 1 



(10) 



Where e^x) = M(x, co t ) - M d (co t ) and {o)i : i = 1, 2, • • • k} is a 
dense set of frequencies which are distributed over in the 
pass-band and stop-band of the filter. A weighting or pen- 
alty error is included to control portions of the actual filter 
response curve that are most important to the filter re- 
sponse. This involves modifying the error to the form: 



ei(x) = W{a) {M{x^i)-MM)) 



(11) 
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Use Rotation matrix to 
generate new points 



Initialization: generate randomly a set of m points 



Figure 4 The spiral optimization algorithm flowchart. 



Evaluation of the fitness values of each point 



Select center to be the point with least fitness 




Term mate 



Where W(co) is a weighting piece wise constant func- 
tion over all frequency space, which is assigned a posi- 
tive value greater than one or less than one to increase 
or decrease the magnitude approximation in a given 
band. 

A recursive filter can be designed by finding a point 
x = x in (11) such that 



e/(#)«0 fori = 1,2, •••/< 



(12) 



Such a point can be obtained when solving the optimi- 
zation problem by minimizing the error function ei(x). 
The design of a recursive filter that approaches a speci- 
fied response M^(c^), can be performed by minimizing 



the error objective function in terms of Lp norm error 
which is formulated as: 



min 



Subject to r a j 
Where 



And 



< 1-6 



CmaxW = max ( e i( x )) 
l<i<l 



(13) 



*»£™{Wf <14) 
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Figure 6 Step response of the filter and time delay (about 100 ms) (SOOF). 



The p is a positive integer. The spiral technique is used 
to minimize fi(x) for increasing values of p. 



In the present design, a constant group delay property is 
considered which is defined as: 



Group-delay objective function 

The group delay is derived from the phase relation, as 
given in equation (4), and is defined as: 



t(x, Co) 



d(p(x 1 co) 
dco 



(15) 



Where (p(x, co) is the phase response of the filter, and 
stated as 



<p(x,co) = arg(H(^ r )) 



(16) 



The group delay for a synthesized recursive filter is de- 
sired to be unchanged in the considered region on one 
hand. On the other hand, a minimum group delay is an- 
other required property needed in certain applications. 



T(pC, 00i 



t c for coi^oor, 



(17) 



Where r c is constant, and co p is the passband regions 
of the filter. In the present application r c is an unknown 
but can be considered as the mean value over the pass- 
band region, which can be determined as: 



— T ( x m, Oi) for 0)iGCO p 
n k ;_i 



(18) 



Where X yyi IS the optimal filter coefficient determined 
by minimizing the magnitude objective function for an 
m-th stability margin parameter. 

The stability margin parameter is varied for discrete 
values, from which an optimal constant group delay is 
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Figure 7 Goup delay response of single objective optimized filter (SOOF). 
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determined by minimizing the following objective error 
function: 

min£ T (m) = ^(r(# w , C0i)-T c ^ m ) for a^Op 

(19) 

The multi-objective optimization problem is solved by 
discretizing the stability margin parameter, the magnitude 
optimization algorithm is used to generate the correspond- 
ing filter coefficients; in which basis the group delay is syn- 
thesized where a feasible and optimal solution can be 
obtained by minimizing the objective function (19). 

The final multiobjective design is obtained by minimiz- 
ing the sum combination of the errors in equations (14) 
and (19). 

The spiral inspired optimization method 

Compared with traditional optimization techniques and 
other global optimizers, the spiral optimization method is 
easy to implement and very efficient in reaching optimum 
solutions. Spiral optimization method has been recently 
developed based on the analogy to spiral phenomena 
(Tamura and Yasuda 2011a; Tamura and Yasuda 2011b). 

Patterns in nature are visible regularities of form found 
in the natural world. These patterns recur in different con- 
texts and can sometimes be modelled mathematically. 



Natural patterns include symmetries, trees, spirals, mean- 
ders, waves, foams, arrays, cracks and stripes. Mathem- 
atics, physics and chemistry can explain patterns in nature 
at different levels. Patterns in living things are explained by 
the biological processes of natural selection and sexual se- 
lection. Studies of pattern formation make use of computer 
models to simulate a wide range of patterns. 

Among the natural patterns, spirals are common in 
plants and in some animals. For example, in the nautilus 
(Figure 2a), each chamber of its shell is an approximate 
copy of the next one, scaled by a constant factor and ar- 
ranged in a logarithmic spiral. 

The spiral phenomena occurring in nature (like the one 
in Figure 2a) are approximated to logarithmic spirals as in 
Figure 2b. Examples of natural spiral dynamics include 
whirling currents, low pressure fonts, nautilus shells and 
arms of spiral galaxies. Logarithmic spirals discrete pro- 
cesses to generate spirals that can form an effective behav- 
iour in metaheuristics. A two-dimensional algorithm has 
been first proposed (Tamura and Yasuda 2011a), and then, 
a more generalized n-dimensional version has been re- 
cently suggested (Tamura and Yasuda 2011b). 

In the present work, the use of the spiral optimization 
technique is presented and used to solve the multiobjective 
IIR filter design. First, based on practical requirements in- 
volved in power system application, a detailed mathemat- 
ical IIR filter design formulation is presented. 



Table 1 Single-objective optimized filter (SOOF) SOS coefficients and gains 



Section Numerator Denominator Gain 

1 1 -1,9520335260 0,9999998000 1 -1,8473812327 0,8836 0,00089772300915 

2 1 -1,9999996000 0,99999980000 1 -1,9999996000 0,999999800000010 1 

3 1 -1,9999996000 0,99999980000 1 -1,8473812327 0,88360000000 1 

4 1 -1,8851094580 0,99999980000 1 -1,8473812327 0,883600000000 1 

5 1 1,99999980000 0,99999980000 1 -1,8473812327 0,88360000000 1 
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Figure 9 Magnitude filter response (MOOF). 



Next, the spiral optimization algorithm detailed and im- 
plemented to solve the optimization task. A first optimi- 
zation example aims at matching a desired magnitude 
response only. Later, the design is improved by considering 
the other requirements including minimum and constant 
group delay. 

Before presenting the n-dimensional spiral optimization 
algorithm, it is worth understanding the two dimensional 
optimization model as some results are just extended over. 

Two-dimensional spiral optimization 

Rotating a point in a 2-dimensional orthogonal coordin- 
ate system (as shown in Figure 3) to the left around the 
origin by 6 can be expressed as: 



: R{6)X 



Where 



Ri{8) 



COS0 

sin(9 



- sin(9 
cos(9 



(20) 



(21) 



Hence, the two dimensional algorithm moves from one 
point to another as: 



xi(k + l) 
x 2 (k + l) 



rR 2 (0) 



%i(k) 
x 2 (k) 



(22) 



Where 0 is the rotation angle around the origin 
(0 < 6 < 2tt) and r is the convergence rate of distance be- 
tween a point and the origin at each iteration k (0 < r < 1). 
The parameter r can be regarded as the "scaling" or radius 
of the logarithmic spiral curvature. 

The spiral model presented earlier has a center only at 
the origin. Hence, it should be extended to have center 
at an arbitrary point x* as: 



x(k + 1) = rR 2 (0)x(k)-(rR 2 (0)-I 2 )x* 



(23) 



Based on the previous formulation, the following opti- 
mization algorithm may be developed: 
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Figure 10 Optimized filter phase response MOOF. 
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Figure 11 Group delay response of optimized filter MOOF. 



♦ Preparation: select the number of search points m > 2, 
the parameters 6 and r and the maximum number of 
iterations k max . 

♦ Initialization: initialize randomly the points; Xi(0) 
i = l..m; in the feasible region and the center x* as 
the point with the least fitness value. 

♦ Updating %{. 



xt(k + 1) = rR 2 (6)x i (k)-(rR 2 (e)-I 2 )x* for i 



l...m. 

(24) 



♦ Updating x*: Select x* as the point with the least 
fitness function in the updated set of points. 

♦ Check for termination criterion: If k = k max then 
stop. Otherwise, start a new iteration. 



n-dimensional spiral optimization 

The extension of the two-dimensional optimization 
algorithm presented earlier is easy to do as one must 



understand how rotation in an n-dimensional space 
is done. Rotation in n-dimension is performed in the 
same way as the two-dimensional rotation taking two 
dimensions at a time. This is defined for dimensions 

U j as: 



COS0 



sin(9 



cos(9 



1 

(25) 



Where the blank elements are zeros. 




0.015 0,02 0.025 
Time (sec) 

Figure 12 Optimized filter step response and time delay (14.9 ms) MOOF. 



u.o4 



Ouadi et al. SpringerPlus 2013, 2:461 
http://www.springerplus.eom/content/2/1/461 



Page 10 of 13 



Table 2 Multi-objective optimized filter (MOOF) SOS coefficients and gains 



Section Numerator Denominator Gain 



1 


1 -1,82933833319 


0,999999800000 1 


-1,0914383222 


0,35402500000 


0,012499922808 


2 


1 -1,4247565639 


0,99999980000 1 


-1,0914383222 


0,35402500000 


1 


3 


1 -1,1469437066 


0,99999980000 1 


-1,1899988100 


0,35402500000 


1 


4 


1 0,227981442276607 


0,999999800000010 1 


-1,09143832226114 


0,354025000000000 


1 


5 


1 -0,467539619942456 


-0,532460226811516 1 


-1,11373644574629 


0,354025000000000 


1 



Hence, there are ^Lii rotation matrices. The resulting 
rotation matrix is then (Tamura and Yasuda 2011b): 

n-l i 

Rn{0) = HOI Rn-i,n+l-i(d)) (26) 
i=l j=l 

Hence the n-dimensional algorithm may be formulated 
similar to the two-dimensional algorithm as: 

♦ Preparation: select the number of search points m > 2, 
the parameters 6 and r and the maximum number of 
iterations k max . 

♦ Initialization: initialize randomly the points; #;(0) 
i= l..m; in the feasible region and the center x* as 
the point with the least fitness value. 

♦ Updating X{. 

Xi(k + 1) = rR n (d)xi(k)-(rR n (d)-I n )x* for i = l...m 

(27) 

♦ Updating: x*: Select x* as the point with the least 
fitness function in the updated set of points. 

♦ Check for termination criterion: If k = k max then 
stop. Otherwise, start a new iteration. 

The flowchart in Figure 4 summarizes the spiral opti- 
mization procedure. 



Results and discussions 

The digital filter to be optimized is to be used to elimin- 
ate harmonics and sub-harmonics in a power network 
with a fundamental frequency of 50 Hz. The filter is a 
bandpass type and is desired to satisfy the magnitude re- 
sponse to ideally pass only frequencies confined in the 
interval [45 Hz, 55 Hz] while rejecting all other fre- 
quency content. The sampling frequency is taken to be 
1800 Hz. The filter is of order 10 and is hence composed 
of 5 cascaded SOSs. At start up, the filter is optimized to 
match the magnitude response specifications only. Next, 
more constraints are added to the optimization process 
including minimum and linear phase and constant group 
delay to enhance the designed filter performance. 

Single objective design 

The purpose of this part is to design a filter which only 
satisfies the magnitude response described earlier with- 
out considering any other performance criteria. The fil- 
ter is designated single objective optimized filter (SOOF) 
hereafter. Figure 5 shows both the desired and the opti- 
mized magnitude response of the digital filter. It can be 
noted that the filter fulfils well the requirements of mag- 
nitude response. Indeed, the filter response falls exactly 
within the desired frequency range [45 Hz, 55 Hz] and it 
attenuates all other frequencies as the overall Sidelobe 
level is lower than -16 dB. 
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Figure 13 Group delay comparison between the MOOF and the classical filter designs. 
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However, the optimized filter is not of practical use as it 
suffers from drawbacks in the dynamic properties. First, as 
shown in Figure 6, the time delay is in the order of 100 ms 
which is not suitable in the present applications as the re- 
quirements specify that the time delay should not exceed 
one cycle (20 ms). In addition, the group delay is not con- 
stant as shown in Figure 7. Furthermore, the phase re- 
sponse is nonlinear as shown in Figure 8 which explains 
the non constant nature of the group delay. As a result, it 
is necessary to include all the preceding performance cri- 
teria into the design process and hence, the problem be- 
comes a multiobjective optimization task. The results are 
summarized in Table 1 where are presented the filter SOS 
coefficients and SOS gains. 

Multiobjective filter design 

The inclusion of the constant and minimum group delay 
in the optimization task besides magnitude response 



criterion produced a filter which satisfies almost all re- 
quirements. The filter is thereafter called multiobjective 
optimized filter (MOOF). The magnitude response of the 
MOOF is shown in Figure 9. The filter magnitude response 
is not as good as the magnitude response characteristic of 
the SOOF. Indeed, the SOOF bandwidth is narrower than 
the MOOF which makes it having better selectivity. 

However, the MOOF performs well in terms of the 
other criteria. In terms of phase response and group 
delay, Figures 10 and 11 illustrate the phase and group 
delay responses of the MOOF. The filter is characterized 
by a linear phase with an almost constant group delay in 
the pass-band contrary to the SOOF phase response and 
group delay illustrated in Figures 7 and 8. Furthermore, 
the group delay is low yielding a small time response as 
shown in Figure 12. Indeed, the time delay of the MOOF 
is about 14.9 ms which conforms to the desired specifica- 
tions. This is due to the fact that the stability margin or 
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equivalently the pole radii have been taken into account 
in the optimization process and these latter have been 
lowered to a value of 0.595 instead of being closer to 
unity. Table 2 summarizes the SOS coefficients and gains 
of the MOOF. 

To validate the obtained results, the MOOF optimized 
filter is compared with the classical filter designs widely 
used in literature. The filters have the same order and 
should meet the same requirements described earlier. 
The filters considered are the least p th optimized filter, 
the elliptic and Chebyshev type II filters. 

The classical filters have magnitude responses that 
meet well the requirements as compared to the MOOF. 
This can be explained by the fact that these filters are 
mainly designed to match a given magnitude response 
which can be considered to be similar to the SOOF opti- 
mized filter. However, according to Figures 11 and 13 
the MOOF has lower and constant group delay over the 
pass-band region compared to the classical filters. In 
fact, the delays of the classical filters are in orders of 200 
and 300 sample delay (the MOOF delay is at 14 samples) 
which is not acceptable for practical requirements. Also, 
the classical design approach relies on matching a re- 
quired magnitude response followed by an all-pass filter 
equalizer. Hence, the classical design adds on the system 
order and the computational complexity. The latter ap- 
proach is not suitable for real time applications. 

To better assess the performance of the filter (MOOF), 
two tests have been performed. In the first, a step sinus- 
oidal signal of 50 Hz is input to the filter and Figure 14 
shows both filter input and output waveforms. It is seen 
that the filter output matches exactly the input except 
for a phase shift and a time delay of less than one cycle. 
Hence, this filter proves to be practical for high speed 
measurement systems where the system accuracy is of 
great importance. In the second test, the previous step 
sinusoidal signal is corrupted with a DC offset, harmonic 
and subharmonic components. The subharmonic com- 
ponent is set to 25 Hz and the harmonics to 100, 150 
and 200 Hz. Both input and output signals are shown in 
Figure 15. As it is clearly seen, the filter succeeded in 
eliminating the DC and harmonic components and miti- 
gating the subharmonic component. 

Conclusion 

The application of the spiral optimization method to de- 
sign a multiobjective digital filter has been considered in 
this paper. The objectives of the filter design were to 
match a desired magnitude response while having a 
minimum and linear phase. At start up, only magnitude 
response has been considered in the optimization task. 
The resulting filter was good in terms of this characteris- 
tic while it showed awful dynamic and phase perform- 
ance. Next, the dynamic properties were included in the 



optimization algorithm to solve a multiobjective task. The 
spiral optimization method has succeeded in attaining the 
optimal design in terms of the previous requirements by 
achieving a compromise between them. The optimized fil- 
ter has been tested and it showed good performance with 
required practical characteristics. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

AO designed the filters using the optimization algorithm, carried out 
research on the filter specifications and characteristics and contributed to 
the theoretical part of the manuscript. HB coordinated the research effort, 
assessed and commented the obtained results and revised the drafted 
manuscript. AR applied the optimization algorithm and searched and wrote 
the rationale behind it. All authors contributed to the final manuscript. All 
authors read and approved the final manuscript. 

Received: 16 June 2013 Accepted: 4 September 2013 
Published: 13 September 2013 

References 

Ababneh Jl, Bataineh MH (2007) Linear phase FIR filter design using particle 

swarm optimization and genetic algorithms. Digit signal process. 

doi:10.1016/j.dsp.2007.05.011 
Antoniou A (1993) Digital filters: analysis and design. McGraw Hill, New York 
Antoniou A (2005) Digital signal processing: signals, systems and filters. McGraw 

Hill-Professional, New York 
Beliczynski B, Kale I, Cain GD (1992) Approximation of FIR and IIR digital filters: an 

algorithm based on balanced model reduction. IEEE Trans on Signal Process 

40(3):532-542 

Chen HC, Chen OT (2006) Particle swarm optimization incorporating a 
preferential velocity-updating mechanism and its application in IIR filter 
design. IEEE Trans Syst Man Cybern:1 190-1 195. SMC'06 October 08-1 1, 2006, 
Taipei, Taiwan 

Cortelazzo G, Lightner MR (1984) Simultaneous design in both magnitude and 

group-delay of IIR and FIR filters based on multiple criterion optimization. 

IEEE Trans Circuits Syst 32(5):949-967 
Dumitrescu B, Niemisto R (2004) Multistage IIR filter design using convex stability 

domains defined by positive realness. IEEE Trans Signal Process 52(4):962-974 
Gordana J-D (2006) Design of IIR notch filters with maximally flat or equiripple 

magnitude characteristics, 14th edition. Europeon Signal Processing 

Conference (EUSIPCO 2006), Florence, Italy 
Ho CYF, Ling BWK, Chi ZW, Shikh-Bahaei M, Liu YQ, Teo KL (2008) Design of near- 

allpass strictly stable minimal-phase real valued rational IIR filters. IEEE Trans 

Circuits Syst II Exp Briefs 55(8):781 -785 
Ikehara M, Funaishi M, Kuroda H (1992) Design of complex all-pass networks 

using Remez algorithm. IEEE Trans Circuits Syst II 39(8):549-556 
Inukai T (1980) A unified approach to optimal recursive digital filter design. 

IEEE Trans Circuits Syst CAS-27(7):646-649 
Jiang A, Kwan HK (2009) IIR digital filter design with new stability constraint 

based on argument principle. IEEE Trans Circuits Syst I Reg Papers 

56(3):583-593 

Jiang A, Kwan HK (2010a) Minimax design of IIR digital filters using iterative 
SOCP. IEEE Trans Circuits Syst I Reg Papers 57(6):1 326-1 337 

Jiang A, Kwan HK (2010b) Minimax design of IIR digital filters using SDP 
relaxation techniques. IEEE Trans Circuits Syst I Reg Papers 57(2):378-390 

Jing Z (1987) A New method for digital allpass filter design. IEEE Trans ASSP 
35(1 1):1557-1564 

Karaboga D (2004a) Simple and global algorithm for engineering problems, 

DE. Turk J Elect Engin 12(1):53-60 
Karaboga N (2004b) Design of minimum phase digital IIR filters by using genetic 

algorithm. Proceedings of the 6th Nordic signal processing symposium- 

NORSIG2004, June 09-1 1, 2004. Espo, Finland 
Karaboga N (2005) Digital filter design using differential evolution algorithm. 

EURASIP J Appl Signal Process 8:1269-1276 
Kit Sang T, Kim-Fung M (1998) Design and optimization of IIR filter structure 

using hierarchical genetic algorithms. IEEE Trans Ind Electron 45(3):481 -487 



Ouadi et ol SpringerPlus 2013, 2:461 
http://www.springerplus.eom/content/2/1/461 



Page 13 of 13 



Krusienski DJ, Jenkins WK (2004) Particle swarm optimization for adaptive IIR filter Tseng CC, Lee SL (2002) Minimax design of stable IIR digital filter with prescribed 

structures. Congr Evol Comput 1:965-970 magnitude and phase responses. IEEE Trans Circuits Syst I Fundam Theory 

Lai XP (2009) Optimal design of nonlinear-phase FIR filters with prescribed phase Appl 49(4):547-551 

error. IEEE Trans Signal Process 57(9)3399-3410 Wang WP, Zhou LF, Qian JX (2004) Fir filter design: frequency sampling filters by 

Lai XP, Lin ZP (2010) Minimax design of IIR digital filters using a sequential particle swarm optimization algorithm. Proc Int Conf Mach Learn Cybern 

constrained least-squares method. IEEE Trans Signal Process 58(7)3901-3906 4:2332-2327 
Lang M, Laakso Tl (1994) Simple and robust method for the design of allpass Wang X, Meng X, He Y (2006) A novel neural networks-based approach for 

filters using least-squares phase error criterion. IEEE Tarns Circuits Syst II designing FIR filters, The sixth world congress on intelligent control and 

41(1):40-48 automation. 1:4029-4032. 1, June 21-23, Dalian, China 

Lang MC (1999) Constrained design of digital filters with arbitrary magnitude and Zhang X (2008) Design of maximally flat IIR filters with flat group delay 

phase responses. Ph.D. Dissertation. Vienna Univ. Technol, Vienna, Austria responses. Signal Process 88:1 792-1 800 

Lang MC (2000) Least-squares design IIR filters with prescribed magnitude and Zhang X, Iwakura H (1999) Design of IIR digital allpass filters based on eigenvalue 

phase response and a pole radius constraint. IEEE Trans Signal Process problem. IEEE Trans Acoustics Speech Signal Process 47(2):554-559 

48(1 1):3109-31 21 

Lertniphonphun W, McClellan JH (2001) Unified design algorithm for complex 

FIR and IIR filters. Proceedings of the IEEE international conference on 

accoustics, speech and signal processing 6. Salt Lake City, Utah, 

pp 3801-3804. May 2001 
Lu WS (1998) Design of recursive digital filters with prescribed stability margin: a 

parameterization approach. IEEE Trans Circuits Syst II Analog Digit Signal 

Process 45(9):1 289-1 298 
Lu WS (2000) Design of stable minimax IIR digital filters using semidefinite 

programming. Proc IEEE Int Symp Circuits Syst 1:355-358. Geneva, Switzerland 
Lu WS (2006) An argument-principle based stability criterion and application to 

the design of IIR digital filters. Proc IEEE Int Symp Circuits Syst:443 1-4434. 

ISCAS2006, Island of Kos: Greece, May 2006 
Lu WS, Hinamoto T (2003) Optimal design of IIR digital filters with robust stability 

using conic quadratic-programming updates. IEEE Trans Signal Process 

51(6):1581-1592 

Lu WS, Pei SC, Tseng CC (1998) A weighted least-squares method for the design 

of stable 1-D and 2-D IIR digital filters. IEEE Trans Signal Process 46(1):1-10 
Lutova A (1997) Design of computationally efficient elliptic IIR filters with a 

reduced number of shift-and-add operations in multipliers. IEEE Trans Signal 

Process 45(10) 

Omoifo Ol, Hinamoto T (2004) Optimal design of stable recursive digital filters 
using unconstrained optimization methods. Proc Int Midwest Symp Circuits 
Syst 2:ll-49-ll-52 

Pan ST (2009) Design of robust D-stable IIR filters using genetic algorithms with 

embedded stability criterion. IEEE Trans Signal Process 57(8):3008-3016 
Peng SC, Chen BS, Chiou BW (1992) Simultaneous design in both magnitude and 

group-delay of IIR and FIR filters based on multiple criterion optimization. IEE 

ProcG 139(5):586-590 
Quelhas MF, Petraglia A (2009) Digital filter design optimization using partial cost 

functions. Proc IEEE Int Symp Circuits Syst:285-288. Taipei, Taiwan 
Recioui A (2012) Sidelobe level reduction in linear array pattern synthesis using 

particle swarm optimization. J Optimiz Theory Appl 153(2):497-512. 

doi:1 0.1 007/sl 0957-01 1-9953-9 
Saab S, Lu WS, Antoniou A (1999) Design and implementation of lowpower 

IIR digital filter systems. Proc IEEE Int Symp Circuits Syst:391 -394. 

Orlando, FL3 

Sanathanan CK, Koerner J (1963) Transfer function synthesis as a ratio of two 

complex polynomials. IEEE Trans Autom Control AC-8(1):56-58 
Sreeram V, Agathoklis P (1992) Design of linear-phase IIR filters via impulse- 
response gramians. IEEE Trans. Signal Process 40(2):389-394 
Storn R (1996) Differential evolution design of an IIR filter. IEEE Int Conf Evol 

Comput:268-273. 20-22 May, 1996, Nagoya, Japan 
Storn R (2005) Designing nonstandard filters with differential evolution. Signal 

processing magazine. IEEE 22:103-106 
Sullivan James L, Adams JW (1998) PCLS IIR digital filters with simultaneous 

frequency response magnitude and group delay specifications. IEEE Trans 

Signal Process 46(11)306-321 
Sullivan JL, Adams JW (1997) PCLS IIR filters with simultaneous frequency 

response magnitude and phase related specifications. Conference record of 

the thirty-first asilomar; conference on signals, systems and Computersl. 

Pacific Grove, California, pp 705-709 
Tamura K, Yasuda K (201 1a) Primary study of spiral dynamics inspired 

optimization. IEEJ Trans electrical Electron Eng 6(S1 ):S98-S1 00 
Tamura K, Yasuda K (201 1 b) Spiral dynamics inspired optimization. J Adv Comput 

Intell Intell Inform 15(8):1 116-1 122 
Tseng CC (2004) Design of stable IIR digital filter based on least P-power error 

criterion. IEEE Trans Circuits Syst I Reg Papers 5 1 (9):1 879-1 888 



doi:1 0.1 1 86/21 93-1 801 -2-461 

Cite this article as: Ouadi et ol.: Optimal multiobjective design of digital 
filters using spiral optimization technique. SpringerPlus 2013 2:461. 



Submit your manuscript to a SpringerOpen 0 
journal and benefit from: 

► Convenient online submission 

► Rigorous peer review 

► Immediate publication on acceptance 

► Open access: articles freely available online 

► High visibility within the field 

► Retaining the copyright to your article 



Submit your next manuscript at ► springeropen.com 



J 



